Effect of Decoherence on the Dynamics of Bose-Einstein Condensates in a Double-well 

Potential 
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We study the dynamics of a Bose-Einstein condensate in a double- well potential in the mean-field 
approximation. Decoherence effects are considered by analyzing the couplings of the condensate to 
environments. Two kinds of coupling are taken into account. With the first kind of coupling domi- 
nated, the decoherence can enhance the self-trapping by increasing the damping of the oscillations 
in the dynamics, while the decoherence from the second kind of condensate-environment coupling 
leads to spoiling of the quantum tunneling and self-trapping. 
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Bose-Einstein condensates(BECs) in a double-well po- 
tential exhibit many fascinating phenomena that are ab- 
sent in thermal atomic ensembles, for example, quantum 
tunneling and self-trapping p], 0, @, 0, H, @]. Quantum 
tunneling through a barrier is a paradigm of quantum 
mechanics and usually takes place on a nanoscopic scale, 
such as in two supperconductors separated by a thin 
insulator [7[ and two reservoirs of superfluid helium con- 
nected by nanoscopic apertures [1,0]. Recently, tunneling 
on a macroscopic scale (jim) in two weakly linked Bose- 
Einstein condensates in a double- well potential has been 
observed [loj. Similar to tunneling oscillations in super- 
conducting and superfluid Josephson junctions, Joseph- 
son oscillations are observed when the initial population 
difference is chosen to be below the critical value. When 
the initial population difference exceeds a critical value, 
an interesting feature of the coherent quantum tunneling 
between the two BECs is observed, i.e., tunneling oscil- 
lations are suppressed due to the nonlinear condensate 
self-interactions. This phenomenon is known as macro- 
scopic quantum self-trapping. 

The interactions between the condensate and noncon- 
densate atoms lead to decoherence. Describing deco- 
herence by fully including the quantum effects requires 
sophisticated theoretical studies that include the effect 
of noncondensate atoms. Treating the noncondensate 
atoms as a Markovian reservoir, master equations that 
govern the dynamics of the condensate atoms might be 
derived [ll], 03, G3 • I n f ac ^ hi the experiments on BECs, 
trapped atoms are evaporatively cooled and they contin- 
uously exchange particles with the noncondensate atoms. 
Thus standard procedure in quantum optics for open sys- 
tems would naturely lead to master equations for treating 
atomic BECs. This Markovian treatment for the BECs 
also can be understood as the presence of lasers for trap- 



ping/detection of atoms, which will polarize the atoms 
and thus couple them to the vacuum modes of the elec- 
tromagnetic field [14]. On the other hand, due to the un- 
avoidable interaction of the BECs with its environment, 
the decoherence is always there in BECs, hence the char- 
acterization of decoherence in this system become inter- 
esting. Because different decoherence may have different 
effects on the dynamics of the BECs, the character of 
decoherence in the BECs may be read out from the dy- 
namics of the BECs. Indeed, as we shall show, different 
BEC-environment coupling leads to different final popu- 
lation imbalance of the BECs in the double-well poten- 
tial. This may be used to characterize the decoherence 
in the double-well systems. 

In this paper, we study the effect of decoherence on the 
dynamics of BEC in a double- well potential by studying 
the evolution of the master equation for the BEC within 
a mean-field framework, where the number of atoms 
in the condensates is supposed to be infinity and the 
quantum fluctuation is negligible. To derive the master 
equation, we need modeling the environment and BEC- 
environment coupling. However, this is not a easy task 
that we do not address at present. Instead, we write 
the master equation by analyzing the effects of environ- 
ment ly induced decoherence. When analyzing decoher- 
ence effects on the dynamics of BECs in a double-well 
potential, we are interested in answering two basic ques- 
tions: (1) What effects are made by the decoherence on 
self-trapping in the BECs? And (2) how does the de- 
coherence affect the quantum tunneling in BECs in a 
double- well potential? 

Consider BECs in a double-well potential, the wave 
function of BECs can be expressed as the supperposition 
of individual wave functions in each well, 



a R \R) + cll\L), 
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where \R) and \L) denote the wavefunction of the right 
and left well, respectively. The coefficients clr and cll 
of the expansion satisfy the two-mode Gross-Pitaevskii 
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equation (GPE) (setting ft = 1), 

i 9(a R \ =H U R \ 
dt\a L J \a L J w 

The Hamiltonian is given by 

/l + |(M 2 -|a L | 2 ) \ \ 

\ \ -l-^W?-\a L ?))' 

(3) 

where 7 is the energy bias between the two wells, c stands 
for the nonlinear parameter describing the condensate 
self-interaction, and V depending on the height of the 
barrier, is the coupling constant between the two conden- 
sates. In this paper, we shall focus on 7 = 0, i.e., the case 
of BECs in a symmetric double- well potential. This situ- 
ation is interesting because the amplitude distributions of 
all eigenstates are symmetric, leading to Josephson oscil- 
lations in the absence of decoherence. With the Markov 
approximation, the master equation that results from the 
condensate-environment coupling takes the form, 

d 

i di P = \- H p>A+ C ^ 

C(p) = i^(2ApA^ - pA^A-A^Ap), (4) 

where T denotes the decoherence rate of the conden- 
sates, and A stands for an operator of the condensates. 
This master equation can be derived by assuming that 
the condensates-environment couplings take the form 
Hj ~ J2j9j(Abj + /i.e.), where gj denotes the constant 
of interaction between the condensates and the environ- 
mental mode by The condensate operator A in general 
is expressed as a linear superposition of three pauli op- 
erators, i.e., A = X x a x + X y a y + X z a z with notations 
a z = \R)(R\ - \L)(L\, a x = \R)(L\ + \L)(R\, and it is 
similar for <j y . The values of X a (a = x,y,z) depends 
on the source of decoherence and its couplings to the 
environment. For example, X y = = is for the en- 
vironment that dephasingly couples to the condensates, 
while X z = is for the environment leading the BECs 
into dissipation. H p takes the same form as in Eq.(j3]), 
except a change \a x \ 2 — > p xx = (x\p\x),x = R,L. 

To start with, we consider the case of A = <j + = 
a x + icjy. This situation happens in the case where the 
double- well potential is formed by using a Raman scheme 
to couple two hyperfine states in a spinor BEC. The con- 
densate in the upper hyperfine states decays into the 
lower one, reminiscent of atomic spontaneous emission. 
The dynamics of the master equation is studied by nu- 
merical simulations, the results are presented in Fig. [1] 
and Fig. [2J In Fig. l-(a) and l-(b), we have plotted the 
population of condensates in the Left well l-(a) and Right 
will l-(b). The initial state is all the condensate atoms 
in the Right well, and the decoherence rate has been set 
to be T = 0.1F. In contrast, the dynamics of the con- 
densate in the Left well without decoherence (T = 0) is 
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FIG. 1: (color online) Populations of the condensates in the 
Left well \L) [(a) and (c)] and Right well \R) [(b)]. The con- 
densates were initially prepared in the Right well, the deco- 
herence rate was chosen r = 0.1V in (a) and (b), while r = 
in (c). A jump- like change at c = 2V in l-(a) and l-(b) clearly 
appears due to the decoherence effect. The nonlinear coupling 
constant was plotted in units of V, while the time in units of 
1/V in all figures in this paper. 
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FIG. 2: The same as FigEQ(b) but with larger r, T = 0.3V 
for the upper panel, while T = 0.5V for the lower panel. 
The jump-like change in the population disappears with T 
increasing. 

presented in Fig. [U-(c). Clearly, the decoherence increase 
the damping of the oscillations. When the nonlinearity 
characterized by c is small compared to the tunneling V, 
the oscillations of the population are suppressed, and the 
condensate finally remains in the two wells with equal 
probability. If the nonlinearity is large with respect to 
the tunneling V and the population imbalance exceeds a 
critical value, the condensate would be locked in one of 
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FIG. 3: Norm of the off-diagonal element of the density matrix 
|pi2 1 as a function of time and c, which is usually used to 
characterize the decoherence. The parameter chosen is F — 
O.iy, and the condensates were initially in the Right well. 



the wells, depending on the initial population. We would 
like to notice that the population change drastically in 
the vicinity of the critical value c = 2V, this is due to 
the suppression of population oscillations by the decoher- 
ence. With decoherence increasing, the jump-like change 
near the critical value in the population becomes unclear, 
as Fig|2] shows. That means the decoherence may deter- 
mine the final population imbalance in the two wells. On 
the other hand, the nonlinear interaction together with 
the initial population and the relative phase can affect the 
decoherence effect, which may be characterized by |pi2|, 
i.e., the norm of the off-diagonal element of the density 
matrix. This effect was shown in figure [3j We would like 
to note that the jump-like change in FigQ] might appear 
at different c, depending on the fixed points around that 
the population imbalance and the relative phase oscil- 
late. For example, our simulations show that the jump- 
like change could appear at c = V with initial relative 
phase = 7r, and non-zero population imbalance @. 

Next, we take A = a x , corresponding to BECs in a 
spatial double-well potential. The tunneling is driven 
by an environment (or by fluctuational fields), leading 
to the decay in the quantum tunneling. An alternative 
BEC system can be formed by using a Raman scheme to 
couple two degenerate hyperfine states in a spinor BEC. 
The driving fields may fluctuate, resulting in decoherence 
in the quantum tunneling. We have performed extensive 
numerical simulations for the master equation Eq. (j4j) . Se- 
lected results, divided into three regimes by the nonlin- 
earity, are presented in Fig. [HE1 and El Figj4] shows the 
dynamics of the condensate in the self-trapping regime. 
The decoherence clearly increases the amplitude of oscil- 
lations in the population first, then increases the damping 
of the amplitude of oscillations, meanwhile it averagely 
decreases the population imbalance, and finally spoils the 
self-trapping. In the self-trapping regime, the frequency 
of the oscillation depends on the nonlinear parameter c, 
the initial population imbalance and relative phase, as 
well as the coupling constant V between the BECs. This 




FIG. 4: (Color online) Population of the condensates changes 
with time. The parameters chosen are, (a) c — 3V, F — 
0.01V, and (b) c = 3V,F = 0. Self-trapping occurs in the 
absence of decoherence, as figure 4-(b) shows. The decoher- 
ence first drives the BECs from the self-trapping regime to 
the quantum tunneling regime, then it destroys the quantum 
tunneling. 
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FIG. 5: The same as FiggJ but with c = V. Solid line is 
plotted for the population of BEC in the Right well, while 
dotted line for the Left. 



can be found by comparing Fig HIE] and El With c and V 
fixed, the decoherence changes the population imbalance, 
this results in the frequency change as shown in Figj4j-(a). 
In the quantum tunneling regime(Fig. [5]), the decoher- 
ence increase the damping of the oscillation as expected. 
And at last, in Fig El we have plotted the dynamics of the 
condensate in the regime between the quantum tunneling 
and self-trapping. We see that the decoherence increases 
the tunneling at the beginning of evolution, and then de- 
stroys the quantum tunneling/self-trapping after a few 
cycle of evolution. 

In summary, we have studied the dynamics of Bose- 
Einstein condensate in a double-well potential. The 
dynamics is govern by the master equations with the 
condensate operator that comes from the condensate- 
environment coupling. Two kinds of decoherence char- 
acterized by cr_|_ and <j x are considered. By numerically 
solving the master equation, we show that there is a 
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FIG. 6: This figure shows the population of the condensate 
at the critical value c = 2V. (a) is plotted for F = 0.01V, 
while (b) is for r = 0. It confirms that the decoherence first 
leads the BEC from the self-trapping regime to the quantum 
tunneling regime, then spoils the quantum tunneling. 



jump-like change in the BEC population due to the first 
kind of decoherence + With the decoherence rate in- 
creasing, the jump-like change in the population becomes 
unclear, resulting in the decoherence-rate dependent self- 
trapping. When the second kind of decoherence (a x ) 
dominated, the decoherence first drives the BEC from the 
self- trapping regime to the quantum- tunneling regime, 
then it destroys the quantum tunneling in the double- 
well system. The limitation of this paper is that we have 
treated the environment Markovian and have ignored the 
quantum fluctuation in the condensate atoms. This may 
limit the application of the formalism to real double- well 
systems. 
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